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The operations drip and mate considered in (mem)brane computing resemble the operations cut and 
recombination well known from DNA computing. We here consider sets of vesicles with multisets 
of objects on their outside membrane interacting by drip and mate in two different setups: in test tube 
systems, the vesicles may pass from one tube to another one provided they fulfill specific constraints; 
in tissue-like P systems, the vesicles are immediately passed to specified cells after having undergone 
a drip or mate operation. In both variants, computational completeness can be obtained, yet with 
different constraints for the drip and mate operations. 



1 Introduction 

One of the basic operations used in the field of DNA computing was introduced by Tom Head in |[20l 
more than twenty years ago, when he formalized the operation of splicing, well-known from biology as 
an operation on DNA strands: given two strings of symbols x and y, the splicing operation consists of 
cutting x and y at certain positions (determined by the splicing rule) and pasting the resulting prefix of 
x together with the suffix of y as well as pasting the resulting prefix of y together with the suffix of x, 
respectively. Formally, if we apply the splicing rule (mi,j<2; 113,114), then the results of splicing x and y are 
z and w where x = x\U\U2X2, y = y\u?,U4y2, and z = X\u\Uny2, w = y\u^U2X2 with Mi ,1*2,1*3, u^,x\,X2,y\ ,y2 
being strings over a given alphabet V. In the case of real DNA sequences, the alphabet consists of four 
letters, i.e., A, C, G, T, representing the four bases adenine, cytosine, guanine and thymine; the cutting is 
realized by restriction enzymes, and the recombination by ligases. 

In ifTTl . the range of Turing machines was encoded using iterated splicing on multisets (sets with 
multiplicities associated to their elements). The splicing operation then mainly was used as a basic tool 
for building a generative mechanism, called a splicing system or H system, as formalized by Gheorghe 
Paun in the following way: given a set of strings (axioms) and a set of splicing rules, the generated 
language consists of the strings obtained in an iterative way by applying the rules to the axioms and/or 
to the strings obtained in preceding splicing steps. If we add the restriction that only strings over a 
designated subset of the alphabet are accepted in the language, we obtain an extended H system. As 
already shown in (H and then in QUI for a class of related systems, in that way we can only obtain 
regular languages. Yet when considering multisets of strings as already done in IfTTl or by adding control 
mechanisms as used in the area of formal language theory (e.g., see [ 10]) like checking for the occurrence 
or the absence of specific subsequences in the strings, then the (extended) H systems were shown to be 
very powerful generative mechanisms, i.e., characterizations of recursively enumerable languages in 
terms of various types of H systems were obtained, for example, see |[25l and lfl5l . 

The idea of computations using test tubes as in HI (Leonard Adleman describes the implementation 
of a small instance of the travelling salesman problem) was formalized to test tube systems using the 
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splicing operation in Q ; again, computational completeness of this computing model could be proved. 

The two subprocesses in splicing, i.e., the cutting by enzymes and the recombination by ligases, were 
introduced as independent operations in cutting and recombination systems; computational completeness 
of several variants of systems using these operations of cutting and recombination instead of splicing was 
exhibited in |[T9lk computational completeness of test tube systems using these operations was proved 
in lfl3lk computational completeness of H systems using cut and paste together with other regulation 
mechanisms as checking for the occurrence of specific symbols or subsequences was shown, too. In 
lfl4l . computational completeness of test tube systems using splicing or cutting and recombination with 
the minimal number of two test tubes was shown; this result is optimal with respect to the number of test 
tubes, because due to Dennis Pixton's results from [30], with one test tube only regular languages can be 
generated. 

For an overview on many interesting models and variants in DNA computing the interested reader is 
referred to the monographs [29] and lf2T6 . 

About ten years ago, another intriguing paradigm based on biology was introduced by Gheorghe Paun 
- membrane systems, soon called P systems (see [26 ]); multisets of objects evolve according to evolution 
rules associated with the membranes arranged in a hierarchical membrane structure. A computation 
consists of transitions from one configuration to the next one, usually applying the rules in a maximally 
parallel manner (i.e., applying a multiset of rules that cannot be extended anymore); the result of a halting 
computation is given by the objects present in the final configuration in a specified output membrane or by 
the objects which leave the external membrane of the system (the skin membrane) during a computation. 
In tissue(-like) P systems (e.g., see |[22l ) the membranes are arranged in an arbitrary graph structure 
instead of a tree structure as in the original model of P systems. A great variety of variants has been 
investigated during the last decade, with the objects being atomic elements or strings, the rules evolving 
these objects and/or moving them through membranes (in P systems) or from one cell to another one (in 
tissue P systems). Many models have turned out to be computationally complete, even with a quite small 
number of membranes or cells, respectively, and with quite restricted variants of rules. The interested 
reader is referred to the monograph [27 ] for an introduction to the wide field of (tissue) P systems and to 
the P systems web page 11241 for the actual state of the art in P systems. 

Whereas in P systems and tissue P systems the objects are placed inside the membranes, in the variant 
of membrane systems introduced by Luca Cardelli (see 0), the objects are placed on the membranes. 
The computations in these models also called brane calculus are based on specific ways to divide and 
fuse membranes and to redistribute the objects on the membranes (e.g., see 0, (3j , (91), the rules usually 
being applied in a sequential way in contrast to the (maximally) parallel way of applying rules in P 
systems. Various attempts have already been made to combine different models from the area of P 
systems and of brane calculi (e.g., see O, (28]). Following this research line by investigating tissue 
P systems with the brane operations mate and drip, in lfl6l computational completeness results were 
obtained both for symbol objects as well as for string objects. As we shall see later in this paper, the 
notations and results given there allow for drawing a close connection to specific models as investigated 
in the area of DNA computing and described above. 

The rest of the paper is organized as follows: After some preliminary definitions, we present our 
definitions for the operations drip and mate and then show the relation of these operations from the area 
of (mem)brane computing with the operations cut and paste used in the area of DNA computing. In 
the fourth and in the fifth section, we prove the computational completeness of test tube systems and of 
tissue-like P systems using drip and mate rules working on sets of multisets. A short summary of results 
concludes the paper. 
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2 Preliminary Definitions 

For the basic elements of formal language theory needed in the following, we refer to any monograph in 
this area, in particular to ||32l . We just list a few notions and notations: N denotes the set of non-negative 
integers (natural numbers), the set of all ^-vectors of natural numbers. By N k RE we denote the set of 
all recursively enumerable sets of vectors of natural numbers. 

V* is the free monoid generated by the alphabet V under the operation of concatenation; its unit 
element is the empty string, denoted by A. The length of a string x € V* is denoted by |jc|; by RE 
(RE (k)) we denote the family of recursively enumerable languages (over a ^-letter alphabet). For any 
family of string languages F, PsF denotes the family of Parikh sets of languages from F and NF the 
family of Parikh sets of languages from F over a one-letter alphabet. In the following, we will not 
distinguish between NRE, which coincides with PsRE(l), and RE (1). 

Let {a i,...,a„} be an arbitrary alphabet; the number of occurrences of a symbol a, in x is denoted 
by \x\ a . \ the Parikh vector associated with x with respect to a\,...,a n is (|jc| , ...,|jc| ) . The Parikh 
image of a language L over {a\,...,a n } is the set of all Parikh vectors of strings in L. For a family of 
languages FL, the family of Parikh images of languages in FL is denoted by PsFL. A (finite) multiset 
(mi,a\) ... (m n ,a n ) with m; 6 N, 1 < i < n, is represented as any string x the Parikh vector of which with 
respect to ai,...,a n is (mi, ...,m n ) . 

In the following we will not distinguish between a vector (mi , . .. , m n ) , its representation by a multiset 
(nt\,a\) ... (m n ,a n ) or its representation by a string x with Parikh vector [\x\ a , \x\ a ) = (mi, ...,m„) . 
In that sense, PsRE (k) = N k RE. 

A deterministic register machine is a construct M = (n,B, Iq, //,,/), where n is the number of registers, 
B is a set of instruction labels, 1$ is the start label, 4 is the halt label (assigned to HALT only), and / is a 
set of instructions of the following forms: 

• Zi : (ADD(r),/2) add 1 to register r, and then go to the instruction labeled by fa; 

• l\ : (SUB(r),Z2,^) if register r is non-empty (non-zero), then subtract 1 from it and go to the 
instruction labeled by 1%, otherwise go to the instruction labeled by ly, 

• lh : HALT the halt instruction. 

A deterministic register machine M accepts a set of (vectors of) natural numbers in the following 
way: start with the instruction labeled by Iq, with the first registers containing the input, as well as 
all other registers being empty, and proceed to apply instructions as indicated by the labels and by the 
contents of the registers. If we reach the HALT instruction, then the input number (vector) is accepted. 
It is known (e.g., see ll23Tl ) that in this way we can accept all recursively enumerable sets of (vectors of) 
natural numbers. In fact, for accepting any L G PsRE (k) we need at most k + 2 registers. 

3 The Operations Mate and Drip 

The reader is supposed to be familiar with basic elements of membrane computing, (e.g., see the mono- 
graph [27 ] and the P systems web page ll24l ). as well as of brane calculi (see, e.g., @). 

The operations we are dealing with in this paper are inspired by the ideas from both areas of P 
systems and of brane calculi: we consider cells with the objects being placed on the membranes of the 
cells (for example, as already considered in lf3TI and [28]) - we will call them vesicles in the following 
- and the operations mate and drip which are taken from the area of brane calculi and very closely 
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related to the model of (mem)brane systems already considered in various papers (e.g., see ||28| , 
[2]), where multisets or strings (in the biological interpretation we may speak of proteins) are placed on 
the membranes. In order to visualize a vesicle with the multiset of objects w assigned to its membrane 
we will use the notation [ ] similar to the notation used in the model of (mem)brane systems. 
The two operations drip and mate we shall use in this paper are defined as follows: 

drip : (u\c\v;y,z) 
mate: (u\a,b\v;x) 

These formal notations describe how to split one cell into two cells (drip) and how to fuse two cells 
into one (mate). 

Following the notations of used in the model of (mem)brane systems these operations are inter- 
preted for the concept of vesicles used in this paper as follows: 

The drip operation (w|c|v;j,z) splits a vesicle (membrane, cell) [ ] sucvw into the two vesicles [ ] 
and [ ] zvw \ (u\a,b\v;x) fuses a vesicle carrying the multiset sua and the vesicle carrying the multiset bvw 
into one vesicle which then has the multiset suxvw, i.e., ab is replaced by x and the remaining multisets 
are taken as they are. In fact, this means that from the two vesicles [ ] and [ ] hvw we get the vesicle 

[ ] suxvw " 

When dealing with strings, the formal notation is exactly the same as given above for the case of 
multisets of objects with the only difference that suy, zvw, and sucvw have to be interpreted as strings in 
exactly the sequence they are written which means that in the case of the drip operation, we start from a 
string sucvw which then is split at the site c yielding the two new strings suy and zvw, hence, s and w are 
not arbitrary anymore. 

In the general case, a,b,c,s,u,v,w,x,y,z can be arbitrary strings over an alphabet V (no matter 
whether these are interpreted as multisets of objects or directly as strings). Computational complete- 
ness for tissue P systems and (mem)brane systems with mate and drip operations working on strings 
using a minimal number of membranes was shown in [18] and H7). 

In contrast to this general case which we shall use in this paper, several restrictions were imposed in 

0: 

1. a,b,c £ V; 

2. b = X,z = X; 

3. V 7^ X, UX ytz X. 

As a special variant of the drip rule dealing with a multiset on the skin membrane of a vesicle we 
also consider the one-sided drip rule where the whole rest of the multiset on the membrane of the vesicle 
to be divided is put to the first target vesicle, i.e., 

dripl : (w|c|v;y,z) 

which in this case means that from a vesicle [ ] sucv we get the two vesicles [ ] and [ ] v ,. 

In contrast to 0, where the weight of a drip rule (w|c|v;y,z) is defined as the length of the multiset 
ucv and the weight of a mate rule (u\a,b\v;x) as the length of the multiset uxv, we here - as already 
considered, for example, in [ 18] - define |«cv;yz| to be the weight of the drip rule (w|c|v;;y,z) and \uabvx\ 
to be the weight of a mate rule (u\a,b\v;x). When using drip rules, one-sided drip rules, and mate rules 
of weight at most k we shall write dripk, dripl^, and mate^, respectively, as parameters in the systems 
(test tube systems and tissue-like P systems) defined in the succeeding sections. 
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3.1 Relating DNA Computing and Membrane Computing 

As already exhibited in Ifl2l . we may observe a coincidence with operations well known from the area 
of DNA computing when looking carefully into the definitions of the operations mate and drip and 
the results of applying them to strings: in |fl9l , the operations cutting and recombination of strings, 
operations which are closely related to the splicing operation, were considered; as we shall exhibit in 
the following lines, cutting respectively its more general variant cut is similar to the operation drip and 
recombination respectively its more general variant paste is similar to the operation mate. 

The cutting operation means cutting a string into two pieces, with adding strings on the cutting sites 
of the cut pieces; the recombination operation means fusing two strings thereby eliminating substrings 
at the fusion sites of both strings. The substrings added at the cutting sites and those eliminated at the 
fusion sites can be interpreted, for example, as electrical charges of molecules. 

More general variants are the cut and paste operations formally to be written as follows: 

cut : («|c|v;y,z) 

cut one string into two strings 
paste : (u\a,b\v;x) 

recombine two strings into one 

The cut operation (w|c|v;j,z) means splitting one string into two strings: a string sucvw is split into the 
two strings suy and zvw, i.e., c is eliminated and replaced by y at the end of the first substring and by z 
at the beginning of the second substring; formally this can be written as sucvw =^ (suy, zvw). The paste 
operation (u\a,b\v;x) means fusing two strings to one string: a string sua and a string bvw are fused 
to the single string suxvw, i.e., ab is replaced by x and the remaining substrings are taken as they are; 
formally this can be written as (sua, bvw) => sucvw. In cutting and recombination systems, we have the 
restrictions x = X and c = A. 

Looking carefully into these notations of the operations cut and paste as well as drip and mate and 
the effect of applying them to strings or multisets, we realize that we have got identical notations: 

mate I paste : (u\a,b\v;x) 
drip/cut : (w|c|v;j,z) 

With respect to the interpretation in tissue P systems with mate and drip operations, a string assigned to 
a cell corresponds with this string itself in the interpretation of DNA computing. Hence, we observe that 
the mate and drip operations and the cut and paste operations are closely related. In that way, results 
established and questions/problems raised for systems using the mate and drip operations may also be 
established/raised for the corresponding systems using the cut and paste (cutting and recombination) 
operations and vice versa. 

As a specific example of relating the two areas of DNA computing and membrane computing, we 
take over the idea of working with sets from DNA computing instead of working with multisets as usually 
done in the area of membrane computing to the model of tissue-like P systems with mate and drip rules. 
On the other hand, we will use the drip and mate rules in test tube systems working on multisets of 
elementary objects placed on membranes. 

4 Test Tube Systems with Drip and Mate Rules 

In this section, we prove our first main result establishing the computational completeness of variants of 
test tube systems with mate and drip rules working on sets of multisets, i.e., as objects in the test tubes 
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we consider sets of vesicles carrying multisets of elementary objects (symbols) on their skin membrane, 
and as operations acting in the test tubes we take the operations drip and mate processing these vesicles. 

We use the following general definition for test tube systems as in lfl4l . where the contents of the 
tubes is redistributed to selected test tubes according to specific filters: 

A test tube system (a TTS for short) a is a construct 

(0,0 T ,n,A,p,D,E) 

where 

1. O is a set of objects; 

2. Oj is a set of terminal objects, Oj C O; 

3. n, n > 1, is the number of test tubes in a; 

4. A = (A\,...,A n ) is a sequence of sets of axioms, where A,- CO,Ki<n; 

5. p is a sequence (pi, ...,p„) of sets of test tube operations, where p,- contains specific operations for 
the test tube 7), 1 < i < n; 

6. D is a (finite) set of prescribed output/input relations between the test tubes in a of the form 
(i,F,j) , where 1 < i < n, 1 < j < n, i ^ j, and F is a (recursive) subset of O; F is called a filter 
between the test tubes 7} and 7); 

7. E C {i | 1 < i < n} specifies the set of output tubes. 

In the interpretation used in [ 14 ], the computations in the system a run as follows: At the beginning 
of each computation step the axioms are distributed over the n test tubes according to A, hence, test tube 
7} starts its first computation step with A,-. Now let L, be the contents of test tube 7]- at the beginning 
of a computation step. Then in each test tube the rules of p, operate on L,, i.e., we obtain p* (Li), 

where p* (Li) = U~ =0 pf , ' > (L,) with p ( - (L,) being defining inductively as follows: pf (L,) = L, and 

pf I+1) (L{) = p\ n) (L{) U Pi (p\ n) (Li) \ for n > 0; for any set L, p,- (L) is the set of all objects obtained by 

applying rules from p,- to objects from L. The next substep is the redistribution of the p* (L,-) over all 
test tubes according to the corresponding output/input relations from D, i.e., if (i,F,j) £ D, then the test 
tube Tj from p ( * (L,) gets p ; * (L,) whereas the rest of p ; * (L,) that cannot be distributed to other test tubes 
remains in Tj. The final result of the computations in o consists of all terminal objects from Oj that can 
be extracted from an output tube f from E , i.e., we take pj- (Lf) n Or- 

In this paper, we allow a more relaxed view of processing the operations in the test tubes and the 
succeeding redistribution of the objects therein, i.e., we assume that at any moment objects fulfilling the 
specific constraints given by a filter (i,F,j) G D may pass from test tube 7} to test tube 7), with some 
copies remaining in 7J. In the limit, the same results can be obtained in that way as in the strict inter- 
pretation as described before, yet our more relaxed interpretation allows for a much easier description of 
development of objects as will be seen in the following. 

The multisets only consisting of terminal objects found on vesicles in an output tube form the set 
of results generated by a test tube system, and the family of all such sets of multisets over a terminal 
alphabet with cardinality k generated by test tube systems using at most m test tubes, axioms of weight 
at most /, drip rules of weight at most q, and mate rules of weight at most p is denoted by 

TTS m (axiomi,drip q ,mate p ) (k) = PsRE (k) . 
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Theorem 1. TTS m (axiomi,mate p ) (k) = PsRE (k) for all m > 3, / > 3, p > 5, k > 1. 

Proof. Let M = (n,B,po,pt,,I)be a register machine with n registers accepting L € PsRE (k) ; moreover, 
let Badd and Bsub denote the sets of labels of the ADD- and SUB -instructions in /, respectively, i.e., 

Badd = {/i|/i:(ADD(r),/ 2 )€/}, 
Bsub = {h \ h '■ (SUB (r) ,h,h) G /} . 

Then we construct a TTS a 

(0,0 T ,3,A,P,D,{3}) 

with three test tubes and mate rules of weight five generating L, with the contents of register i represented 
as the number of symbols b { as follows: 

The objects in O are vesicles of the form [ ] w with w being a multiset over an alphabet V to be 
specified below; yet we may simply represent such an object by any string representing w; hence, we can 
also write O = V* where 

V = B U {X,Y,Z,F} U {a t \ \ <i<k}yj{bi\ \ <i<n) 
U {A t \ I € Badd} U {A u A'iA" I / G Bsj/b} • 

In the same sense, we will write Oj = Vj with Vj = {a, \ 1 < i < k}. 

In the first test tube T\ , we initialize the simulation of a computation in the register machine M with 
obtaining (vesicles carrying) multisets of the form Xa" 1 ...d^b" 1 ...b n £ using the axioms 

{X,Zlo}U{aibiY | 1 < i<k} 

and the mate rules (X \,Y |;) and (X |,Z | lo;); with applying the second rule, we start the simulation of 
a computation in the register machine M. 

Moreover, l\ : (ADD(r) , h) £ I is simulated by the axiom A; t hb r and the mate rule (X 1 1\ ,A^ \ ljb r \ ) . 

For/i : (SUB(r),/2,/3) G /, the subtract case is simulated using the axiom A/, h and the mate rule (X | 
l\b r , A/, |/2;)- The case when we guess the contents of register r to be zero is started with using the 
axiom A J together with the mate rule (X \ A^;). The computation is then continued in test tube T2 
where the rule (X \ Aj f ,A'[ \h;) with the axiom A" I3 allows for sending back the multiset in case that the 
guess has been correct. Appearance checking (testing that no symbol b r is present) in the zero case for 
l\ : (SUB(r),/2,^) G / is accomplished by the corresponding filter in 

(l,Ui< r <„(V r U{X}U{bi I \ <i<n,i^r}U {A' h \ h : (SUB(r),/ 2 ,/ 3 ) € /})*,2) 

from test tube T\ to test tube J2 and 

(2,(V -{A' u A'l\l^B SUB })\\) 

from test tube T2 back to test tube T\ . 

The terminal results are collected in test tube T3 by eliminating the symbol X which is present in every 
multiset representing a configuration of a computation in M as soon as the final label If, has appeared with 
using the mate rule (| hX,F |; ) with the axiom F in test tube T\ and then letting these terminal multisets 
get through the filter (!,{«, | 1 < i < k}* ,3) from test tube T\ to test tube J3. 
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The sets of axioms, rules, and prescribed output/input relations (filters) A, p, and D, respectively, can 
easily be collected from the descriptions given above: 

A = (A b A 2 ,0), 

Ai = {X,Zh,F}U{aibiY\l<i<k} 

U {A h l2br\h:{ADD(r),l 2 )€l} 

U {A h l 2 ,A' h \h : (SUB(r),/ 2 ,/ 3 ) €/}, 

A 2 = {A'lh \h : (SDB(r),/ 2 ,/ 3 ) €/}, 

P = (Pl,p2,0), 

Pi = {(X|,y |;),(X|,Z|/o;),(|/*X,F|;)} 

U {(X I Z 1 ,A /l |/ 2 fe r ;)|/i:(ADD(r),/ 2 )€/} 

U {(X | /^ r , A,, | / 2 ;),(X | Z x ,| A' h ;) \ h : (SUB(r),/ 2 ,/ 3 ) € /} , 

p 2 = {{X\A' h ,A>[ |/ 3 ;)|/i:(SUB(r),/ 2 ,/ 3 )€/}, 

D = {(l,Ui< r <„(y r U{X}U{^- | 1 < / < tvV r} U {A^ |/i : (SUB(r),Z 2 ,Z 3 ) €/})*, 2), 
(2, (V - {A^Af I I G , 1), (1, {a,- | 1 < i < k}* ,3)}. 

As desired, we use only three test tubes, axioms of weight at most three, and mate rules of weight at 
most five; moreover, the filters in the prescribed output/input relations of the TTS a are of the very special 
and simple form (i,W*,j) with W C V or finite unions of such filters. These observations complete the 
proof. □ 

As an alternative to having all the axioms in the test tubes as indicated in the proof constructed above, 
we may use the single axiom g and the drip rule 

(UI;A,) 

for each axiom A. Hence, we immediately obtain the following result: 

Corollary 2. TTS m (axiomi,drip q ,mate p ) (k) = PsRE (k) for all m > 3, / > 1, p > 5, q > 4, k > 1. 

Proof. All required axioms can be computed from the single axiom g by using the drip rule (| g |;A, ) - as 
well as by using (| g \;g,) for g itself - in each of the two test tubes T\ and T 2 . As a small technical detail 
we mention that the computations in these new test tube systems need an additional step at the beginning 
to initialize the two test tubes T\ and T 2 with the corrsponding set of axioms. □ 

Another interesting variant is the use of one-sided drip rules instead of mate rules: looking carefully 
into the proof of Theorem 1 and the mate rules used there we realize that the second vesicle always 
carries an axiom. In general, if bv is the whole second vesicle, then the mate rule 

(u | a, b | v;x) 

can be simulated by the one-sided drip rule 

(u | a | ; vx, ) , 

i.e., we put everything to the first vesicle and thus in fact obtain only one result by applying this rule. 

Corollary 3. TTS m (axiomi,drip\ q ) (k) = PsRE (k) for all m > 3, / > 1, q > 4, k > 1. 

Proof. According to the proof of Corollary 2, we can get every axiom by a one-sided drip rule. Moreover, 
as explained above, every mate rule (u\a,b\ v;x) used in the proof of Theorem 1 can be replaced by the 
one-sided drip rule (u \ a \;vx, ). □ 
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5 Tissue-like P Systems with Mate and Drip Rules 

In this section, we prove our main result establishing the computational completeness of variants of 
tissue-like P systems with mate and drip rules working on sets of multisets. 

A tissue-like P systems with mate and drip rules (tP system for short) IT is a construct 

(y,V T ,n,A,R,i ) 

where 

1 . V is a finite set of symbols; 

2. Vj is a set of terminal symbols, Vt C V; 

3. n, n > 1, is the number of cells in IT; 

4. A = (Ai, ...,A„) is a sequence of sets of axioms, where A; C V*, 1 <i <n, describing the initial 
contents of the cells; 

5. R is a set of rules of the form 

Tr.r^Tj 

with i,j G {/ | 1 < / < n}, i ^ j, and r being a drip or mate rule over V; 

6. j'o G {/ | 1 < / < n) specifies the output cell. 

A computation in IT starts with the initial configuration described by A; a computation step then 
consists of applying the rules 7} : r — > Tj in the i-th cell - the application of a rule 7] : r — > 7) means 
applying r to objects in (the source) cell 7} and sending the resulting vesicle(s) to (the target) cell 7) - in 
a maximal way in that sense that every vesicle that can undergo the application of a rule will be affected 
by a suitable rule, yet as we are dealing with sets of vesicles, this also means that any vesicle or any pair 
of vesicles has to be used with every possible rule by which it can be affected. 

The multisets only consisting of terminal objects found on vesicles in the output cell i'o form the 
set of results generated by IT, and the family of all such sets of multisets over a terminal alphabet with 
cardinality k generated by tissue-like P systems using at most m cells, axioms of weight at most /, drip 
rules of weight at most q, and mate rules of weight at most p is denoted by 

tP m (axiomi , drip q , mate p ) (k) = PsRE (k) . 

Theorem 4. tP m (axiomi, drip q , mat e p ) (k) = PsRE (k) for ailm>5,/>3,/?>5,g>5,£>l. 

Proof. Let M = (n,B,po,Ph,I) be a register machine with n registers accepting L € PsRE (k); then we 
construct a tissue-like P system IT 

(y,V T ,S,A,R,S) 

generating L. We start with the following initial vesicles in the five cells: 

A x = {B s \ s e {X,Zl ,F} U {aibiY \ 1 < i < k} 
U{A h l 2 b r \h:(kDD(r),l 2 )eI} 
U{A h k,A' h |/i:(SUB(r),/ 2 ,/ 3 )G7}}, 

A 2 = ©, 

A 3 = {E r l 3 , F r D r | h : (SUB {r),l 2 ,h) €/}, 
A 4 = {A r \h:{Sm(r),l 2 ,h)eI}, 
A 5 = 0. 
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In general, for generating a multiset s in the first cell 7\ we use the following rules in T\ and T 2 : 

7i : (| B s Ij^B^Bjj) -»• T 2 , 

T 2 : (| R,B' S B S |; ) — > Ti generates s in T\, 

T2 : (| R,B' s s |; ) — > Zi regains fi s in 7\. 

Moreover, for sending back from T2 to T\ a multiset containing the specific symbol X indicating a 
multiset on a vesicle representing a configuration of a computation in the register machine M, we use the 
special symbol R with the rule 

T 2 :(X\,R\;)^T h 

For the initialization as already explained in the proof of Theorem 1, we take s = X, s = aibjY for 
1 < i < k, and s = Zl$ as well as the rules 

7i:(X|,y |;)->r 2 and 

7i:(X|,Z|Jo;)->r 2 ; 

with applying the second rule, we start the simulation of a computation in the register machine M. 
For simulating an ADD-instruction l\ : (ADD(r),/2) E / we take s = AiJ 2 b r and the rule 
T l :{X\l u A h \hb r ;)^T 2 . 

For simulating a SUB -instruction l\ : (SUB(r),/2,^) El in the case that subtraction is possible we 
take s = A/j h and the rule 

T\:{X\ hb r ,A h \ h\)^T 2 . 

In all the cases described so far, the main work is done by a rule of the form T\ : r — > T 2 using a rule 
in T\ with the result being sent to cell T 2 , where with the application of the rule 

T 2 :(X\,R\;)^Ti 
we already described before, the result is sent back to cell T\ . 

For simulating a SUB -instruction l\ : (SUB(r),/2,^) E I in the case that subtraction is not possible 
we take s = A\ y and guess that no b r occurs, but now send the result to cell T3 : 

T l :{X\h,\A' h ;)^Ti. 

Checking for the occurrence of b r now is accomplished by the following rules affecting a vesicle 
containing X within a cycle of 2; in even computation steps, the rule T3 : (| B r ,X \ b r \) — > T4 "kills" 
vesicles containing b r by sending them to cell T4 thereby also erasing the symbol X so that it cannot be 
affected by a rule anymore. If no b r occurs, then one step later the rule T3 : (| A J ,E r \ Z3; ) —> T 2 sends the 
vesicle with the desired label I3 back to cell T\ via cell T 2 (hence, in total the simulation of this case takes 
four steps). The symbols A r ,B r ,C r and D r ,E r ,F r , respectively, allow for having the desired checking 
symbols B r and E r in T3 at the right moment, i.e., if a vesicle has "survived" B r , then E r will finish the 
simulation of the zero-case of the SUB-instruction. 

T 4 :(\A r \;B n C r A r )^T 3 , 

T 3 : (\B n C r |A r ;) ^ T 4 , 

T 3 :(\B n X\b r ;)^T 4 , 

T 4 : (| D r \;E r h,F r D r ) T 3 , for h ■ (SUB(r) ,l 2 ,h) G h 
T 3 :(\E r l 3 ,F r \D r ;)^T A , 

T 3 : (\ A' h ,E r \ l 3 ;)^T 2 , for h:(SUB(r),l 2 ,l 3 ) EI. 
To obtain the output vesicles in T5, we apply the rule 

7i:(/AX|,F|;)->r 5 . 



Rudolf Freund & Marian Kogler 



133 



In sum, we obtain the following set of rules R: 

R = {T 1 :(\B s \;R,B , s B s s)^T 2 , 
n:(\B s \;R,B' s B s s)^T 2 , 

T 2 :(\R,B' s s\;)^n 
| s e {X,Zl ,F} U {mbiY | 1 < i < k} 
U{A h l 2 b r \h:{km(r),h)eI} 
U{A h l 2 ,A' h \h : (SUB(r),Z 2 ,Z 3 ) £/}} 
U {n : (X \,Y |;) -> T 2 ,T X : (X \,Z \ Z ;) -> T 2 } 
U {r 2 : |;)->7i,7i : (Z A X |,F |;)->r 5 } 

U {Ti : (X | l u A h | Z 2 Z? r ;) -> T 2 | h : (ADD(r),/ 2 ) G /} 
U {r 1 :(X|/ 1 ,|A; i ;)^r 3 ,r 1 :(X|/ 1 fc r ,A, 1 |/ 2 ;)->T 2 , 

r 3 : (| |fe r; ) -> r 4 ,r 3 : (| A^£ r I h;) -> r 2 , 

T 4 : (I A r |;S r ,C r A r ) -> T 3 ,r 3 : (| B r ,C r |A r ;) -> T 4 , 
T 4 : (I D r \;E r l 3 ,F r D r ) -> T 3 ,r 3 : (| £ r / 3 ,F r | D r ;) -> T 4 , 
T 3 : (|A; i5 £ r | Z 3 ;) -> T 2 \ h : (SUB(r),Z 2 ,Z 3 ) € /} 

We emphasize once more that the simulation of any computation step of the register machine M takes 
an even number of steps (i.e., two or four), and also in the initial phase, i.e., the generation of the axioms 
and the initial configurations Xwlo with w G {afti \ 1 < i < k}* in the first cell T\ takes an even number 
of steps, which guarantees that the zero-check performed by the interplay of rules in the cells T 3 and 
T 4 works correctly. Finally, we mention the computation in IT never stops and every element of L will 
appear as the multiset on a vesicle in the output cell at some moment during the computation in IT and 
will be sent to cell T5 again in each odd step of the computation after its first appearance in T$, as every 
computation of the register machine M can be started again after any even number of computation steps 
in n. These observations complete the proof. □ 

6 Conclusion 

As in DNA computing, we have considered sets of objects instead of multisets as mostly considered in 
the area of P systems. The operations cut and recombination well known from DNA computing have their 
counterparts as the operations drip and mate considered in (mem)brane computing. We have investigated 
the computational power of specific variants of the operations drip and mate on sets of vesicles with 
multisets of objects on their outside membrane acting in test tube systems, where the vesicles pass from 
one tube to another one provided they fulfill specific constraints, and in tissue-like P systems, where the 
vesicles are passed to specified cells after having undergone a drip or mate operation. In both setups, we 
have proved computational completeness, even with different variants of the drip and mate operations. 
As far as the descriptional complexity of the test tube systems with respect to the number of test tubes 
and of the tissue-like P systems with respect to the number of cells and in both cases with respect to 
the weight of the mate and drip operations is concerned, improving the obtained results in these respects 
remains as a challenging task for future research. 

Acknowledgements: The first author gratefully acknowledges many interesting discussions with Gheo- 
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